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Abstract. Using dynamical systems theory and a detailed numerical analysis, the late-time behaviour of 
tilting perfect fluid Bianchi models of types IV and Vllf, are investigated. In particular, vacuum plane-wave 
spacetimes are studied and the important result that the only future attracting equilibrium points for non- 
inflationary fluids are the plane-wave solutions in Bianchi type VII^ models is discussed. A tiny region of 
parameter space (the loophole) in the Bianchi type IV model is shown to contain a closed orbit which is 
found to act as an attractor (the Mussel attractor). Prom an extensive numerical analysis it is found that 
at late times the normalised energy-density tends to zero and the normalised variables 'freeze' into their 
asymptotic values. A detailed numerical analysis of the type VII^j models then shows that there is an open 
set of parameter space in which solution curves approach a compact surface that is topologically a torus. 



1. Introduction 

In recent years there has been much analysis of spatially homogeneous (SH) perfect fluid cosmological 
models with equation of state p = (7—1)/!, where 7 is constant Pj. For the class of tilted (non-orthogonal) SH 
models , the Einstein field equations have been written as an autonomous system of differential equations 
in a number of different ways ^ ^ Tilted Bianchi models contain additional degrees of freedom, and 
new dynamical features emerge with the inclusion of tilt. A subclass of models of Bianchi type V 
and models of Bianchi type II TU' have been studied. 

In a local stability analysis of tilting Bianchi class B models and the type VIq model of class A were 
studied. The asymptotic late-time behaviour of general tilted Bianchi type VIq universes was analysed in 
detail in |12| . and this was generalized to Bianchi type VIq cosmological models containing two tilted 7-law 
perfect fluids in Finally, the tilting Bianchi type VI/j and Vllh models were recently studied in 

(self-similar tilting solutions ^2] of Bianchi type Ylh were also studied in [TB]'). 

In particular, the late-time behaviour of Bianchi models of type IV and VII/i was analysed in jl4| . with an 
emphasis on the stability properties of the vacuum plane- wave spacetimes. It was shown that for 2/3 < 7 < 2 
there will always be future stable plane-wave solutions in the set of type IV and VII^ tilted Bianchi models 
and it was proven that the only future attracting equilibrium points for non-inflationary fluids (7 > 2/3) are 
the plane-wave solutions in Bianchi type VII/i models. A tiny region of the parameter space was discovered 
in the type IV model in which a closed orbit occurs. Moreover, for the type VII^j models a Hopf-bifurcation 
occurs and more attracting closed orbits appear. 

In this paper we will use a dynamical systems approach to analyse the tilting Bianchi type IV and VII/j 
models in more detail, with an emphasis on numerical analysis. Since this current paper is a companion 
paper to 15 , we shall not reproduce all of the details of the analysis and we shall frequently refer to equations 
in the original paper. However, we shall attempt to make the current paper as self-contained as possible. 

We note that the Bianchi Vllh models are of special importance in cosmology. We recall that Bianchi type 
Vllft, models are of maximal generality in the space of all SH models and are sufficiently general to account 
for many interesting phenomena. In particular, these cosmological models are the most general models that 
contain the open Friedmann-Robertson- Walker models and are consequently of special interest. Therefore, 
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these anisotropic generalizations of the open Friedmann-Robcrtson- Walker models have been studied in 
detail in order to examine the dynamical behaviour of exact solutions of the Einstein field equation and, in 
particular, to study the evolution of large scale anisotropics in models that contain the standard low density 
universe as a limiting case. 

There are a number of observational constraints on Bianchi models. For example, the structure of the 
Cosmic Microwave Background radiation anisotropy over intermediate scales can be used to constrain the 
open Bianchi type VII/i cosmological models ^12, . If the universe is open, then a distorted quadrupole or 
hotspot effect can occur jj^. However, the resulting constraint is of no practical use if the density lies close 
to the critical density (as suggested by recent observations). However, even if the density lies close to the 
critical density, in Bianchi VII/i models the anisotropic 3-curvature produces a general relativistic geodesic 
spiralling effect The spiral plus quadrupole are both focussed by the hotspot effect allowing a possible 
discrimination of the open Bianchi VII/i models. 

The asymptotic and intermediate evolution of non-tilting perfect fluid Bianchi VH/j models has been 
studied using a detailed qualitative analysis 1 , in conjunction with numerical experimentation jl8j . It is 
known that although Bianchi VII/i universes can isotropize at late times, the set of models that do so are 
of measure zero in the the class of all Bianchi models |2ijJ. The occurrence of intermediate isotropization 
in this class of non-tilting universes has been studied in 21 using approximate solutions of the Einstein 
field equations. More recently the evolution of the shear in Bianchi VII/i models was discussed heuristically 
|22j : in particular, it was claimed that if initial conditions at the Planck time consistent with the Planck 
equipartition principle are imposed, then the shear will decay sufficiently by the present epoch in order to 
be consistent with the low anisotropy in the Cosmic Microwave Background radiation. In 18 it was found 
that while intermediate isotropization may occur, it does not necessarily do so, and the approximate and 
heuristic analyses of [2] and [55] were critically examined. The linear stability of Bianchi type VII/i vacuum 
plane- wave solutions has recently been studied in 123] . 

The paper is organised as follows. First, we shall write down the equations of motion using the orthonormal 
frame formalism in the W-gauge'. Then, in section|31 we present all of the equilibrium points and investigate 
their stability properties. In section 01 we prove the existence of closed loops and provide some criteria for 
their stability. Then we will study the Bianchi type VII^ (and IV) models numerically in detail. Some of this 
numerical analysis is distributed throughout the paper; however, some explicit examples (e.g., the physically 
interesting case of dust and radiation) arc given in Appendix^ In particular, using a detailed numerical 
analysis of the type VII/i models, we shall fully justify the claim in that there is an open set of parameter 
space in which solution curves approach a compact surface that is topologically a torus. 

2. Equations of motion 

The companion paper |14j contains all the details regarding the determination of the evolution equations 
for the tilted cosmological models under consideration. These equations, written in gauge invariant form, 
allow one to choose the gauge that is best suited to the application at hand. Here, we shall adopt the 
iV-gauge in which the function Nx is purely imaginary; this is ensured by the choice (p' = ^/SXT,^, where A 
is defined by iV = Alm(Nx). The evolution equation for N can then be replaced with an evolution equation 
for A, which ensures a closed system of equations. We note that this choice of gauge is well suited for the 
study of the dynamics close to the vacuum plane- wave solutions ^Ij. For a qualitative analysis of these 
models, the iV-gauge is preferable since the resulting dynamical system is well defined in the Bianchi Vlh 
case. However, care must be exercised when analyzing the Bianchi VII^ case, in which the dynamical system 
is not well defined over the entire region of state space, in particular the variable A could diverge. 

In the notation of |14j . the expansion- normalised anisotropy and curvature variables used in this paper 
are 



Si 1112 + 'iSl3, Sx=S-+iS23, Nx=iN. 
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We will also adopt the dimensionless time parameter r, which is related to the cosmological time t via 
dt/dr = (l/H), where H is the Hubble scalar. 

Using expansion-normalised variables, the equations of motion in the A'^-gauge are (see |14| for the complete 
derivation of the equations): 

(2.1) s; = {q-2)j:+ + 3{j:l, + j:l,)-2N^ + ^{-2vl + vl + vl) 

(2.2) !]'_ ^ (q-2-2V3I]23A)S_ + \/3(E?2-Si3) + 2A7V+^|^^(v|-f2) 

2tr + 

(2.3) S'i2 = (g-2-3S+-\/3S_)Si2-\/3(E23 + S-A)Si3 + ^^|^wi«2 

(2.4) E'i3 = (g-2-3S+ + \/3S_)Si3- V3(E23-S_A)Si2 + ^^wiW3 

(2.5) S'23 = (g - 2)S23 - 2%/3A^2^ + 2V3AE^ + 2V3E12E13 + ^^^^f2W3 

(2.6) N' = + 2S+ + 2V3E23A) iV 

(2.7) A' = 2V3S23(1-A2) 

(2.8) A' = ((7 + 2S+)A 

The equations for the fluids are: 

^[2q - (37 - 2) + 2jAvi + [2g(7 - 1) - (2 - 7) - 7^] V^} 
{T + 2S+) vi - 2%/3Si3W3 - 2%/3I]i2^;2 -A(vI + vl) - %/37V (vl - vl) 



(2.9) 


n' 


(2.10) 




(2.11) 


v'2 


(2.12) 




(2.13) 


v 



VIV2 



T-Y.+ - VSE- j V2-Vi (S23 + 5]_ A) 1^3 + V3AiVwiV3 + (.4 + VsiVj 

T - E+ + VSE-) W3 - \/3 (E23 - S_ A) 1^2 - y/l\NviV2 + VSiv) W1W3 

F(l - V^) 



l-(7-l)y2 



[(87 - 4) - 2(7 - l)Avi - S] 



where 



^ 2 i + (7-i)y2 

Y? = + E^ + E?^ + S^g + 
S - E,bc"c^ c"ca = l, f" = Fc", 

[(37 - 4) - 2(7 - l)Avi] (1 ^ V^) + (2 - 7)^25 
l-(7-l)y2 

G+ = l + (7-l)l^'. 



T 
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These variables are subject to the constraints 

(2.14) 1 

(2.15) 



2Ei A + 2E_7V 



(2.16) 

(2.17) 
(2.18) 







I]i2(iV + \/3A) + Ei3AiV 
Si3(iV- V3yl) + Ei2AiV 



G+ 



The generaUzed Friedmann equation, 



The parameter 7 wiU be assumed to be in the interval 7 € (0, 2) 
(|2.14|l . yields an expression which effectively defines the energy density O. Therefore, the state vector can 
thus be considered X = E_, E12, S13, E23, A, A, vi,V2, W3] modulo the constraint equations (equations 
(|2.15l2.18|l 'l. Thus the dimension of the physical state space is seven (for a given value of the parameter h). 
Additional details are presented in |14| . 

The dynamical system is invariant under the following discrete symmetries ; 



4>i 



[E+,E_,Si2,Ei3,E23,A^, A, A, t!i,'(;2,i^3] 1-^ 
[i;+,i;_,Si2,Ei3,S23,A^, A, A,z;i,t;2,?;3] 1-^ 
[i;+,S_,i;i2,Si3,S23,A^, A, A,i;i,W2,W3] ^ 
[i;+,I]_,i;i2,Si3,S23,A^, A, A,wi,i;2,W3] ^ 



S-, i;i2, Si3, S23, -N, A, "A, -vi,-V2, -vz] 

Si3, S12, E23, -N, A, A, vi,Vi,V2\ 
S-, tSi2, ±Ei3, -E23, N, -A, A, wi, =FW2, iwa] 
!]_., -E12, -Ei3, S23, A^, A, A, wi, -U2, 



These discrete symmetries imply that without loss of generality we can restrict the variables A > 0, A > 0, 
and A'^ > 0, since the dynamics in the other regions can be obtained by simply applying one or more of the 
maps above. The fourth symmetry listed implies that one can add one additional constraint on one of the 
variables Si2,Si3,U2 or W3; however, in general there is no natural way to restrict any one of the variables, 
and hence we will not do so here. Note that any equilibrium point in the region z)2 > has a matching 
equilibrium point in the region V2 < 0. 

2.1. Invariant sets. In this analysis we will be concerned with the following invariant sets: 

(1) T{VIIh): The general tilted type Vllh model. Given by 1 < |A|. 

(2) r(/y): The general tilted type IV model. Given by A = ±1, A > 0. 

(3) T(V): The general tilted type V model. Given by A = 0, A > 0, A = ±1. 

(4) Till): The general type II model. Given by A = ±1, A = 0. 

(5) B{I) Type 1: N = A = V = 0. 

(6) Ti{A): One-component tilted model. Given by W2 = W3 = 0. 

(7) dT{I) "Tilted" vacuum type I: f7 = A = A = 0. 

We note that the closure of the sets T{IV) and T{VIIh) are given by 



Tiyilh) = T{VIIh)UT{II)UB{I)UdT{I), 
(2.19) T{IV) = T{IV)UT{V)UT{II)U B{I)UdT{I). 

Since the boundaries may play an important role in the evolution of the dynamical system we must consider 
all of the sets in the decomposition (|2.19|) . 

3. Future Asymptotic Behaviour 

3.1. Qualitative Analysis. In order to study the future asymptotic behaviour of tilted type IV and VII/i 
models, as a first step we will consider all the equilibrium points of the closed sets T{IV) and T{VIIh)- 
Fortunately, much work has already been completed in this regard. For instance, the Bianchi tilted type II 
models were studied in some detail by Hewitt, Bridson and Wainwright ^O!- They found that the future 
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asymptotic state is either a flat Friedman-Lemaitre model (0 < 7 < 2/3) , a non-tilted Bianchi type II model 
(2/3 < 7 < 10/7), an intermediately tilted Bianchi type II model (10/7 < 7 < 14/9), a line of tilted points 
for 7 — 14/9 or an extremely tilted Bianchi type II model when 14/9 < 7 < 2. Similarly, in a paper by 
Hewitt and Wainwright |5] in which the irrotational Bianchi type V models were analysed, it was found 
again that the flat Friedman-Lemaitre model is stable for 0<7<2/3,a non-tilted Milne model is stable 
for 2/3 < 7 < 4/3 and an intermediately tilted Milne model is stable for 4/3 < 7 < 2. Therefore, if we are 
interested in determining the future asymptotic state of the general tilting Bianchi type IV and YUh models, 
we need only determine whether the corresponding equilibrium points are stable in the full Bianchi IV or 
Vll/i phase space. In particular, all equilibrium points in dT{I) are Kasner circles and related equilibrium 
points. None of these equilibrium points are stable into the future and consequently we will not list these 
below. On the other hand, the set dT{I) is essential in the analysis of the past asymptotic behaviour. We 
should also point out that for the equilibrium points with A'^ = a different gauge has been used in order to 
determine their stability. 

3.1.1. B{I): Equilibrium points of Bianchi type I. 

(1) T{T): [0, 0, 0, 0, 0, 0, A, 0, 0, 0, 0] where A is an unphysical parameter. This represents the flat Friedman- 
Lemaitre model. 
Eigenvalues: 

The remaining equilibrium points are all in dT[I). 

3.1.2. T{II): Equilibrium points of Bianchi type II. All the tilted equilibrium points come in pairs. These 
represent identical solutions (they differ by a frame rotation); however, since their embeddings in the full 
state space are inequivalent, two of their eigenvalues are different. All equilibrium points have an unstable 
direction with eigenvalue —2^/3T,23 corresponding to the variable A. 

(1) CS{II): [^,0,0,0, ^^"^^;!'^ ,^ ^^~g^^^^~^\ 1,0, 0,0,0], 2/3 < 7 < 2. This is the Collins-Stewart 
perfect fluid Bianchi II solution. 
Eigenvalues: 




3(37-2) 3,, ^J^^ ^-^^-^)!^;-^M,-^[x2], ^%-^»)[x2] 



(2) n,{II): [iz^,-i^y/ ai7-io)(2-.)(7,-io) ^o^o^_^ 

7 < 2. This is Hewitt's tilted Bianchi II model with = (37-4)(97^-207+i2) ^ ^ ^ m^2l(^ljli^±fm ^ - 



(37-4)(77-10) „^ J ^2 _ 1 
(ll7-10)(57-4) '^^^'^ e - i. 

Three of the eigenvalues are M22zl}^^ M2^2l^ 3(77-10) ^ .^-^^^i the remaining four eigenvalues being 
the roots of a nasty equation. 



(3) n,{II): [27^,0,|a,fa,-f (57-6),f ^(iH)^l^,l,0,0,6^,-e^], 10/7 < 7 < 2. 
This is Hewitt's tilted Bianchi II model where S^, fl and V are the same as for 0.2(11), and 



_ _ VI . / (ll7-10)(2-7)(77-10) j^jjjj g2 ^ ^ 



4 Y I77-I8 

Three of the eigenvalues are M^t^lLH ^ li52z^^ _ 3(77^40) ^ with the remaining four being the roots 
of the same nasty equation as for 7ii(//). 



(4) Ad/): hi - V5n27/.-+2)(a4- si^^ _ Vi^(^T.-+4)(i7-27,^) _ ^ ^ ^^/g^ 



3 ' 342 7 2^' 2^' 9 

„„„ n ^ ^;2 ^ _ . . . . 

11 ^ ^ (17-27/3^)(64-81/3^) 



This is an intermediately tilted line bifurcation where < < = n 9 7 kY- s ^"^^ 
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n = 013(1458/3^ - 1215/32 + 472). 
Three of the Eigcnvahics are 0, |, |. 



(5) A(77): [0, - i(a-/3), yi W+4)(rr-2r,^) ^ sVM Z, _^] ^ ^ ^ 

14/9. This is an intermediately tilted line bifurcation where < /S^ < a = = 
V27/?2 + 2, ^2 = fflg^gSg) , = 614(1458/34 - 1215/3^ + 472) and = 1. 
Three of the Eigenvalues are 0, |, — |. 

(6) £1(1/): [-i,-^^^,^,-^,-^,^,l,0,e,0,0],0<7<2. This is an extremely tilted Bianchi 
II, with O = 20/57, = 28/57 and = 5^ = 1. 

Eigenvalues: 

~3(2^-~7r^ ' ^±^\/-20615±456x/3n. 

(7) f^a/): [0, (1 + ^) , (1 - ^) , 1' 0> -^'5^. ^'^^l^ < 
7 < 2. This is an extremely tilted Bianchi II, with Q = 20/57, = 28/57 and e'^ = 6'^ = 1. 



Eigenvalues: 



-2(97-14) -11 



, 4/3, -2/3, — ± —y -20615 ± 456V311. 



3(2-7) ' ' ' 3 57 



3.1.3. T{V): Equilibrium points of Bianchi V. 

(1) I{V): [0,0,0,0,0,0,1, A, 0,0,0], < A < 1, 7 = 2/3. Represents a FRW model with a 7 = 2/3 
fluid. 

Eigenvalues: 

0[x2], -2[x5]. 

(2) MiV): [0,0,0,0,0,0,1,1,0,0,0], < 7 < 2. This represents the Milne universe. (S^ = 0,f2 = 

0, = 0) 

Eigenvalues: 

0, -2[x2] ,-(37-2), (37-4)[x3]. 

The zero eigenvalue arises due to the fact that this is part of the one-parameter family of plane- wave 
solutions. 

(3) M{V): [0, 0, 0, 0, 0, 0, 1, 1, 5^^, 0, 0], 6/5 < 7 < 2. This represents an "intermediately tilted" Milne 
model (E2=0,f^ = 0,y2 = gl^). 

Eigenvalues: 

^ ^2(7-1)^ ^' (97-10)(7~l) 
Again the zero eigenvalue is associated with there being a non-isolated line of equilibria. 

(4) M±iV): [0,0,0,0,0,0, 1,1, ±1,0,0], < 7 < 2. This represents "extremely tilted" Milne models 
(S2 = 0,^^ = 0,1/2 = 1). 

Eigenvalues: 

M+ : 0[x2], -2[x2], l[x2], 2, : 0, -2[x2], -4, ^i^I^, _i[x2]. 

2-7 

(5) n{V): [S+, 0,0, 0,0, 0,1,(1 + 1]+), 1,0,0], -1 < E+ < 0, < 7 < 2. This represents a type V model 
with a null fluid. 



FUTURE ASYMPTOTIC BEHAVIOUR OF TILTED BIANCHI MODELS OF TYPE IV AND VII^ 



7 



Eigenvalues: 

0[x2], -2(l + S+)[x2],(l-2S+)[x2], 2(1 - 2S+). 

The two zero eigenvalues are associated with this being part of a two-parameter family of equilibrium 

points {TZ{IV) see below). 

3.1.4. T{IV): Equilibrium points of Bianchi type IV. 

(1) £{IV) : [E+, V-I]+(l + S+), 0, 0, 0, v'-S+(l + S+), 1, (1 + 0, 0, 0], -1< S+ < 0, < 7 < 2. 
These represent 'non-tilted' vacuum plane waves. 

Eigenvalues: 

(1 + S+) ± 2i\/3iv] , -4E+ - (37 - 2), (37 - 4) + 2S+, (37 - 4) - S+ 

(2) £±(/y):[S+, V-S+(l + S+),0, 0,0, ^-^+(1 + ^+), 1, (1 + S+), ±1, 0, 0], -1< E+ < 0, < 7 < 2. 
These represent 'extremely tilted' vacuum plane waves. 

Eigenvalues: 



£+ : 0, -2 
£_ : 0, -2 



(1 + E+) ± 2iV^N] , 0, 2(1 - 2E+), (1 - 2E+)[x2]. 

(l + S+)±2ix/3ivl, -4(1 + S+), -^i^I_^±M±),_(i + 4S+)[x2]. 
J 2-7 



(3) CilV) : [S+, V-S+(l + E+), 0, 0, 0, V-S+(l + 1, (1 + ^^2^±g^,0, 0], -1< S+ < 0, 
6/(5 + 2S+) < 7 < 2. These represent 'intermediately tilted' vacuum plane waves. 
Eigenvalues: 

2-7, 



0, -2 



(1 -I- S+) ± 2iV3N 



7-1 



(1 - 2S+), 



(1 - 2E+)(57 - 6 + 27S+)(37 - 4 + 2S+) (1 - 2S+)(37 - 4) 



x2]. 



(^_1)(9^_10 + 4S+(1-E+)) ' 2(7-1) 
(4) ^±{IV): v/-S+(l + S+), 0, 0, 0, \/-S+(l + S+), 1, (1 + S+) Izll^ 



where ?;2 = ± 



^2(37-4)(-l+2S+)(37-4-S+) 
2(3-27)(l+S+) 



(3-27-l)(S++l)' 

1 < S+ < 0, (4 S+)/3 < 7 < min[4/3, 3/(2 - S+)]. 



These represent 'intermediately tilted' vacuum plane waves. 
Eigenvalues: 

(l-2S+)(57-6) 



0, 



(1 + S+) ± 2iV3N 



3-27 



, 0, Ae, A7, 



where 



AeAr = — 
Ag + A7 = 



(1 - 2S+)2(4 - 37)(37 - 4 - S+)(3 - 27 + 7E+) 
(3-27)G(7,S+) 
(1-2E+)F(7,S+) 
4(3 - 27) (1772 - 4O7 + 24)G(7, S+) 



and 



(3.1) 



G(7,5]+) ^ (57 - 6)S2+ - (18 - 257 + 97')S+ - 3 + 27, 

i^(7,S+) = [2(177^-407 + 24)2+ -337^ + 1217^- 1527 + 66]^ 

-9(7 - if (1217"* - 7367-'' + 16647^ - 16567 + 612) 

Here is G(7, S+) < in the whole region under consideration. F(7, S_|_) = 0, defines a line from 
(4/3,0) to (4/3,-1/4). 
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(5) £±{IV): ^-^+{1 + 0, 0, 0, + 1, (1 + 1^2, -V2I 



where V2 = ^ 



V2(4£++1)(2S+-1) 



vacuum plane waves. 
Eigenvalues: 



6S4 



-1 < < —1/4, < 7 < 2. These represent 'extremely tilted' 



0, -2 



(1 + S+) ± 2is/lN 



(1-2S+)(1 + 2S+) _(1-2S+)(1+4E+) 2(l-2E+)(27-3-7S+) 



S+ ' ' 3S+ ' 3S+(2-7) 

(6) n{IV): [S+, (l-£)^-S+(l + S+),0,0,0, + 1, (1 + S+), 1,0,0], -1< S+ < 0, 

0<£<1, 0<7<2. These represent non-vacuum plane-waves with a null fluid. 
Eigenvalues: Same as for C+{IV). 

3.1.5. T{VIIh): Equilibrium points of Bianchi type Vllh- For all of the vacuum equilibrium points the group 
parameter, h, is given by 

3/iS+(l- A') = (1 + S+), 

and we have defined h = l/Vh. 

(1) CiVIh) : v/-S+(l + 0,0,0, ^-^+{1 + A, (1 + S+), 0,0,0], -1< E+ < 0, 1< A, 
< 7 < 2. These represent 'non-tilted' vacuum plane waves. 
Eigenvalues: 

11+) ± 2iV3N 



0, 



(1 



-4S+ - (37 - 2), (37 - 4) + 2E+, (37 - 4) - S+ ± i/i(l + 



(2) £±(F//^) : V-S+(l + S+), 0, 0, 0, V-E+(l + E+), A, (1+S+), ±1, 0, 0], -1< S+ < 0, 1< A, 
< 7 < 2. These represent 'extremely tilted' vacuum plane waves. 
Eigenvalues: 



£+ : 0, -2 [(1 -I- S+) ± 2z\/3ivl , 0, 2(1 - 2S+), (l-2S+)[x2]. 



(1 + S+) ± 2iV3N] , -4(1 + S+) 



2(57 - 6 + 27S+) 

2-7 



-(1-|-4E+)±2i/i(1 + Eh 



(3) C{VIh) : V-S+(l + S+), 0, 0, 0, + A, (1 + ^^^I^^^g^, 0, 0], -1< E+ < 

0, 1 < A, 6/(5 + 2E+) < 7 < 2. These represent 'intermediately tilted' vacuum plane waves. 
Eigenvalues: 

2-7, 



(1 -I- E+) ± 2i\/3N 



7-1 



(1 - 2E+)(57 - 6 + 27E+)(37 - 4 + 2E+) (1 - 2E+)(37 - 4) 



±ih{i-vi){i + i:+). 



(^_l)(9^_10 + 4S+(l-S+)) 2(7-1) 

(4) n{VIh): [S+,(l-^)V-S+(l + S+),0,0,0,(l-^)V-S+(l + E+),A,(l + S+),l,0,0],-l<E+< 
0, < £ < 1, 1 < A, < 7 < 2. Non-vacuum plane-wave with a null fluid. The group parameter is 
given by 3/iI]+(l - A2)(1 - ^) ^ {1 + 
Eigenvalues: Same as for C+{VIIh)- 

3.2. Local stability of equilibrium points. We quickly observe that none of the equilibria in the Bianchi 

type I and Bianchi type II invariant sets are late-time attractors for the Bianchi IV or Yllh models when 
2/3 < 7 < 2. For the Bianchi type V model, only the isotropic Milne universes can act as future attractors. 
However, as can be seen, these equilibrium points are the isotropic limits of the plane-wave equilibrium 
points of Bianchi type IV and can be extracted directly from this analysis. Furthermore, the isotropic limit 
of the type Vllh model can also be directly extracted from the type VII?^ analysis (even though A diverges 
here). 
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Figure 1. Stability of the Plane- Wave Vacuum Solutions. The regions where the different 
plane- wave equilibrium points are attractors for the Bianchi type IV with V2 + V3 > 0. In the 
left figure, all the boundaries along the left edge marks the instability of the energy density, 
J7. The right figure is a magnified region of the left one showing the loophole (figure taken 
from ^^). 




-0.5 




-0.25 -0.2 -0.15 -0.1 -0.05 



The stability analysis of the type IV model can be summarized as follows: For 2/3 < 7 < 2, the only 
equilibrium points that are future attractors in T{IV) are: 

(1) CilV): stable for < 7 < ^±§±. 

(2) C-{IV) : stable < 7 < 2, -i < E+ < 0. 

(3) T±{IV) : T-{IV) stable in {v2 + I'a) > ^ for max (|, < 7 < min (70, 

(4) E±{IV) : £-iIV) stable for - i < S+ < -i, 5^ < 7 < 2. 

Here, we have defined 70 for any given value of E+ as F{ja, E+) = 0, where F{'j, E+) is defined in ea. H3.1|l . 
Note that there is a region where there are two co-existing future attractors; namely, the region where 
< 7 < 7o. Also, there is a tiny region 70 < 7 < 5_|_2s+ (from now on called the 'loophole') which does 
not contain any stable equilibrium points (see Figure^. 

The stability analysis of the equilibrium points of type VII^ can be summarized as follows: For 2/3 < 
7 < 2, the only equilibrium points that are future attractors in T{VIIh) are: 

(1) C{VIIh): stable for < 7 < 

(2) C-{VIIh) : stable for ^:^:§^ < 7 < 2, -i < 1!+ < 0. 

^JF_|_(7V') is stable in the same region of 7 for the half {v2 -I- fa) < 0. 
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Figure 2 . The figure below displays the dynamical behaviour for values of the parameters 
near the loophole for Bianchi type IV non- vacuum tilted models. The value of 7 is 1.325 
and the initial conditions are chosen to be close to the plane wave solutions for illustrative 
purposes only. A cross indicates the locations of £^{IV), diamonds indicate the position of 
C-(IV) while an asterisk indicates the initial condition. The equilibrium point J^^{IV) is 
located in the center of the spiral and the center of the closed loop. Note how the dynamical 
behaviour depends on the final value of S]-|_. 




3.3. The Bianchi type IV loophole. As pointed out in ^l], in the limiting vacuum case (i.e., limn^o), the 
stability of the plane wave attractors changed as a function of both the parameter 7 and the terminal value 
of E+, henceforth denoted as E^j.. As we have seen from the above there also exists a region of parameter 
space (here we consider the terminal value T,'^ as a parameter) which does not contain any stable plane- 
wave equilibrium points (the loophole). With the use of the bifurcation diagram and through subsequent 
numerical experimentation, a limit cycle was found for values of 7 and E!j_ inside the loophole, see Figure 
n This limit cycle acts as an attractor and will be referred to as the Mussel attractor (see section 5 and 
Figure 2 in m|). Figures [3 and Ol indicate how the dynamical behaviour of the Mussel attractor changes for 
7 = 1.325 and different limiting values of E_|_. 

We note that the limit cycles observed in Figure [21 above are actually a continuum of limit cycles side 
by side. The limiting value of S+ is the determining factor in the final asymptotic state of the system. Of 
course one might argue that the likelihood of actually falling into this family of limit cycles is small. Again 
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Figure 3. The figure below displays the dynamical behaviour for values of the parameters 
near the loophole for Bianchi type IV non-vacuum tilted models. A cross indicates the 
locations of £-{IV), diamonds indicate the position of £-{IV) while an asterisk indicates 
the initial condition. The equilibrium point J-'^(IV) is located in the center of the spiral and 
the center of the closed loop. Note how the dynamical behaviour depends on the value of 
Each trajectory represents the vacuum limit of the orbits observed in the figure above. 
The circular arc represents the boundary of the phase space. 
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through some numerical experimentation it can be shown that there is a set of Bianchi IV initial conditions 
far from the plane-wave solutions of non-zero measure that have these limit cycles as future attractors for 
values of 7 inside the loophole (see Fig. 1 151 where this fact is illustrated for the Bianchi type VII^ loophole). 

3.4. Behaviour of the System in General. One of the most important points and argument in this and 
the companion paper is that at late times f2 ^ and the variables Sab, N, etc., 'freeze' into a particular 
asymptotic value. We have not been able to prove this rigorously, however, there is much evidence to support 
the following: 

Conjecture 3.1. Except for a set of measure zero, non-inflationary (2/3 < 7 < 2j tilted perfect fluid Bianchi 
type IV and Vllh models will at late times approach a vacuum state; i. e. 

lim n^Q. 



There are several results that support this. First, a local stability analysis shows that there are always 
stable vacuum solutions. So the conjecture is at least valid locally. Second, we know that there are no 
non-vacuum equilibrium points acting as local attractors; the preceding analysis excluded this. Third, in 
our fairly extensive numerical analysis no evidence for any other behaviour has been found. In particular. 
Figure 121 shows how the variable E+ freezes in to a particular value of S+. In Appendix IXI more examples 
of this behaviour are shown and we can see explicitly how — > and S-). — > S^j.. 

The local analysis also shows that the variable Sab, N, etc. 'freeze in' to particular values S*^, N* . After 
this freezing has occurred, the system of equations effectively decouples and the equations governing the tilt 
velocities separate. In effect, the tilt equations reduce to the much simpler form Vi = F{vj; S^, SI, ■ ■ • ), 
and hence, we can, at late times, treat all the variables, except vi, V2 and V3, as constants. This 'freezing 
in' process makes these models amenable to analytical investigations, which was done in the analysis 
provides numerical justification for this freezing process (e.g., see Figure|21l. In fact, the numerical investi- 
gations show that the time scale of this freezing process is typically much shorter than many of the other 
features present; e.g. the oscillation in the type IV loophole (see Appendix IXI) . In the Appendix we have 
included a sample of numerical integrations. For the plots we have chosen the particularly interesting values 
7=1 (dust), 1.325 and 4/3 (radiation). In all of the plots we can see how quickly the energy-density, Jl, 
tends to zero and the variable S+ freezes in. On the other hand, the velocities typically take longer to settle 
into a particular asymptotic state (as explicit examples of this, see Figures ITTl and IT^ . We should point out 
that the numerical integrations are done in the full 7-dimensional state space; however, for some figures the 
reduced system is used for illustrative purposes only. 

4. Dealing with Closed Loops 

The dynamical system contains closed periodic orbits. Some of these closed curves act as limit cycles and 
take the role as attractors for the dynamical system. However, closed loops are difhcult to find and describe 
in full generality. 

In our case, the 'freezing in' process enables us to consider the reduced system of equations describing the 
behaviour of the tilt velocities in a vacuum plane- wave background, see |14| . We utilize the identity 

and define {x, p, 6) by 



(4.1) 



X = vi, p = V2+vl, 21)2^3 = p sin 6*, (^2 ~ ^3) = PC0s6'. 
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Then we have for the reduced system 

x' = (T + 2I];)x - {A* + %/3iV* cos B)p, 
p' = 2{T-Y.\+A*x-acose)p, 

(4.2) 9' 2a(A* +sine'), 

where a ~ V3(l ^ x)N* and, by use of the discrete symmetry ^4, we can assume that Q < 9 < 27r. 
Furthermore, these variables are bounded by 

(4.3) 0<p, V'^ = x^ + p<l. 

An asterisk has been added to the variables to emphasize that these should be thought of as the limit values 
for the full system. 

We will utilize the following fact: Assume that c(r) is a (solution) curve, then if / is a function in the 
state space, 

(4.4) r f{ciT))dT^ficit2))-f{c{h)). 

Jtl 

In particular, if c(t) is a closed periodic orbit, 

(4.5) j f\c{T))dT = Q. 

For a periodic orbit, c, with period we introduce the average of a variable B: 

(4.6) {B) = ^<f>BdT. 

We can also say something about the stability of a closed periodic orbit. For example, consider the evolution 
equation for fi, which we write as 17' = rJAn. Assume that c(t) is a closed periodic orbit with period T„. 
Then for every tq and e > 0, there exists a solution curve, c(r) = c(r) + 5{t), for which 

(4.7) \5{T)\<e, To<T<To+Tn. 

(This follows from Proposition 4.2, page 104, in [H].) Hence, the curve c(r) can be used to approximate the 
closed curve c(r). Using c(t) we get 

n{To + T„ 



(4.8) In 



dr^ / Aodr = r„(An)+0(e). 

J To 



n 



Since e can be arbitrarily small, we can approximate 

(4.9) n{To + T„) « n{To) exp(r„ (Ao)). 

Hence, if (Ao) < 0, then ri(ro + Tn) < ^{tq) for a sufficiently small perturbation; i.e., the closed curve is 
stable with respect to il. If, on the other hand, (Ao) > then the closed curve is unstable with respect to a 
perturbation of fl. This result is the corresponding local stability criterion for periodic orbits. 

By calculating (1 — a;) (in two different ways) we are able to find an expression for r„ for orbits in the 
type Vllh model: 

(4.10) r„ = 



(i + s;)(i-(x))' 



where n is an integer. This expression involves the group parameter h, and we note that a larger h implies 
longer orbital period. 
We also need, 
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Lemma 4.1. Assume that there is a closed properly^ periodic orbit, c{t), for the dynamical system \4-S^ - 
Then 

Proof. For (A*)^ > 1, 6* is a strictly monotonic function, so in this case we can use 6* as a variable along the 
curve. Furthermore, 0{tq + T„) = 9q + 2?™, where n G Z, which gives 

2a cos 9dT = / — — -de = In (A* + sin 61) 1^+^™ = 0. 

Since a — ^/3N*{1 — x), the Lemma now follows for (A*)^ > 1. For A* = 1 the equation for 9 implies that 
the periodic curve must be in the invariant set sm9 = —1. Thus the Lemma is trivially satisfied in that 
case. □ 

Theorem 4.2. Assume that there is a closed properly periodic orbit, c{t), for the dynamical system \4--'/^ - 
Then either 

^A^'^^ r l\ (87 - 4) - S; (1 - 2S;)(57 - 6) 

^"^^' (l + I]^)(3-27) ' ^^''^^ 3^^; ' 

or V — 1 and 

1 + , , (1-2S1)(1 + 2S1) 

(4.14) £: (^>-^^. {^n)='- ^. 

Proof. Assume first that there exists a closed periodic orbit with y < 1 and p > 0. Then, using the evolution 
equation for V, we get 

(4.15) (37 - 4) ifiV)) - 2(7 - 1)A* (xfiV)) - {fiV)S) = 0, 

for any analytic function /(V") on (0, 1). Using Lemma [4. II the equation for p' yields 

(T) - s; + A* (x) = 0. 

Using eq. H4.15|l with f{V) = V'^/G- we obtain (T) = (S). Using eq. 14.15|l with f{V) = 1 we can solve for 
(a;), which yields 

^_ (37-4)-!]; 

^""^ (i + s;)(3-27)- 

A similar calculation yields the desired expression for (Ao) . 

For V ^ 1 (which is an invariant subspace) , we have p — 1 — a;^ , and we can rewrite the equation for p: 

ji-x^ y 

1 — a;^ 

Dividing by 2a; and integrating yields (x) — A* as desired. The expression for (Ao) is found analogously. 

Now consider the invariant subspace p = 0. This subspace is homeomorphic to the closed interval [—1,1], 
and hence, no proper periodic orbit can exist there. 

The Theorem now follows. □ 

This theorem states exactly when a closed periodic orbit, if it exists, is stable with respect to the variable 
il. Equation H4.9|l implies that small deviations from = will effectively decay exponentially if (Ao) < 0. 

So far we have only given necessary criteria for closed periodic orbits. However, in some cases we can 
even prove existence of such orbits. 



(4.16) 2^ = 2a; {-'SxY.*^ + A*+2a cos ( 



^By properly periodic we will mean a periodic orbit which is not an equilibrium point. Equilibrium points are trivially 
periodic in the sense that c(r) = c{t + Tn) for any T„ . 
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Theorem 4.3. For every S!j_ G (—1,-1/4), (A*)^ > 1, and 7 G (0,2) there exists one, and only one up to 
winding number, closed periodic orbit, £{VIIh), with V = 1 for the dynamical system i4.S^ . 

Proof. Existence: Given G (—1, —1/4), > 1, and 7 G (0, 2), we consider the invariant subspace V = 1. 
This extremely tilted set is topologically a sphere, S^. For (A*)^ > 1, the variable 9 is strictly monotonic 
and can thus be considered as a "time" variable. On the sphere S'^, there are two, and only two, equilibrium 
points; the North pole, x — 1 and the South pole, x — —1. For S^J_ G (—1,-1/4), both of these equilibrium 
points have two unstable directions, and hence there exist open neighbourhoods Un and Us around the 
North and South Pole, respectively, such that S'^ — {Un ^ Us) is a future trapping region. The compact set 

— {Un U Us) is diffeomorphic to a compact annulus in R^, which allows us to use the Poincare-Bendixson 
theorem. This says that there must be either a future attracting equilibrium point, a future attracting closed 
orbit, or an arbitrary union of points and curves. However, we know there are no other equilibrium points; 
i.e., there must be a closed curve. Existence is thus proven. 

Uniqueness: Since 6 is strictly monotonic, these periodic curves have to be winding around the a;-axis. 
Assume there exists two such curves, c{9), and c{9). Since these two curves cannot cross each other, \x{9) — 
x{9)\ > 0, and equality holds if and only if c = c (up to winding number). By considering the x-coordinate 
for these two closed period orbits, we obtain after a straight-forward manipulation, 

\{x)~{x)\>0, 

where equality holds if and only if c = c. Using Theorem 14 . 21 we have (x) = (2:); hence, c = c (up to winding 
number) and uniqueness is proved. □ 

Theorem 4.4 (Mussel attractor). For A* = 1 and every and 7 taking values in the type IV loophole ( 
7o < 7 < 5^2S* ) ^^c^fi exists a closed periodic orbit, C{IV), for the dynamical system \J~^ 

Proof. In this case we can consider the 2-dimensional invariant subspace sin^ = —1. Again we can use 
Poincare-Bendixson theory, however, in a slightly more complicated form. By a consideration of the equilib- 
rium points and the separatrices the existence of a closed curve can be shown (see | 25| for details and exact 
formulation of the theory). □ 

4.1. Stability of the Bianchi type VII/i plane waves. We now return to the general Bianchi type VII^i 
models. In addition to the analytical results stated above, we have done a comprehensive numerical analysis 
in the full 7-dimensional state space which is summarised in the Figures to follow. The above results, and 
from the comprehensive numerical analysis, indicate that the nature of the late-time attractors, as we change 
the parameter A* from A* = 1 to A* > 1, change according to 



IV 


VII,, 




£{IV) H- 


C{VIIh) 


(equilibrium point) 


£_(/!/) H- 


C{VIh) 


(equilibrium point) 


^-(iv) ^ 


Hvih) 


(closed orbit) 


£^{IV) H- 


- £{vih) 


(closed orbit) 


C{IV) H- 


- T{VIh) 


(torus) 



These closed curves, and the torus, are illustrated in Figures 0] and Analysis indicate that the stability of 
these attractors is preserved under the transition IV 1 — > VII^i. For the curve £{VIIh), this can be shown 
using similar techniques as above; the region of stabilty in terms of the limiting value and 7, is exactly 
the same as for £-{IV) (see Figure P). 
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However, the torus case has to be treated with care; the late-time behaviour is described by a nowhere 
vanishing flow on a torus. The future behaviour of such integral curves come in two classes according to the 
nature of the future limit sets 

(1) Rational curves: The attractor is a closed periodic curve which can be described by its homotopy 
class, 7J, which takes values in the fundamental group of the torus, 7ri(T^) = Z x Z. 

(2) Irrational curves: The attractor is an everywhere dense curve on the torus. 

The rational curves immediately lend themselves to the analysis above; if TJ = (m, n), where n is defined 
as the winding number of 9, then T„ — , where (x) is given by Theorem 14. 21 For an irrational 

curve, c(t), we can utilize the following fact. For any tq and S > 0, there exists a closed curve c(t) with 
period T„ (not an integral curve) such that |c(t) — c(t)| < S/2 for tq < t < tq + T„. In particular, this 
means that |c(ro) — c{tq + T„)| < S, which implies that the limiting curve c(r) can be arbitrary close to a 
closed curve. This implies that we can get arbitrary close to the values given in Theorem 14.21 and hence, we 
get similar restrictions on the limit cycle. The stability analysis for the variable f2 can then be applied in a 
similar manner. Note that for most parameter values in the type VII^ loophole we expect these irrational 
curves to be attractors. Since the flow on the torus is described continuously by the parameters 7 and S^, 
we would expect that locally the direction of the flow takes all real values in a certain interval. Since the 
rational numbers only form a set of measure zero we expect that "most" curves will be irrational. 

5. Conclusion 

We have used a dynamical systems approach and a detailed numerical analysis to analyse the late-time 
behaviour of tilting Bianchi models of types IV and VII/j. We established the equations of motion in the 
'TV-gauge' and we found all of the equilibrium points and investigated their stability, with an emphasis on 
the vacuum plane- wave spacetimes. In particular, it was shown that for 2/3 < 7 < 2 there will always be 
future stable plane-wave solutions in the set of type IV and VII^ tilted Bianchi models and it was proven 
that the only future attracting equilibrium points for non-inflationary fluids (7 > 2/3) are the plane- wave 
solutions in Bianchi type VII^ models. 

The future asymptotic behaviour of tilted 7-law perfect fluid models of Bianchi type IV and Vllh can be 
summarised as follows'^: 
Bianchi type IV:. 

(1) < 7 < 2/3: Asymptotically FRW. 

(2) 2/3 < 7 < 6/5: Asymptotically a vacuum plane wave. Tilt velocity tends to zero. 

(3) 6/5 < 7 < 4/3: Asymptotically a vacuum plane wave. Tilt velocity is asymptotically zero, interme- 
diately tilted, intermediately tilted and oscillatory (in the loophole), or extreme. 

(4) 4/3 < 7 < 2: Asymptotically a vacuum plane wave. Tilt velocity is asymptotically extreme. 

Bianchi type Vllh- 

(1) < 7 < 2/3: Asymptotically FRW. 

(2) 2/3 < 7 < 6/5: Asymptotically a vacuum plane wave. Tilt velocity tends to zero. 

(3) 6/5 < 7 < 4/3: Asymptotically a vacuum plane wave. Tilt velocity is asymptotically zero, inter- 
mediately tilted and oscillatory (either J- (Vllh) or the torus attractor), extreme, or extreme and 
oscillatory. 

(4) 4/3 < 7 < 2: Asymptotically a vacuum plane wave. Tilt velocity is asymptotically extreme or 
extreme and oscillatory. 

The stability of the plane- wave solutions is a key result. 



■^Models with < 7 < 2/3 are subject to the 'no- hair' theorem (originally due to Wald for a cosmological constant I2(j| '). 
The tilted version was proven in |14|. 
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Figure 4. The figure below displays the dynamical behaviour for values of the parameters 
near the loophole 7 = 1.325 for Bianchi type Vllh non- vacuum tilted models (for the figures 
h — 1 is used). When = —.2637, the future asymptotic behaviour is dominated by 
the stable closed curve £(VIIh)- As increases to a value of — —.2578, the closed 
loop £{VIIh) becomes unstable (still exists), and a new stable closed loop appears in the 
form of T{VIIh)- When Y,*j^ — —.2518, the future asymptotic behaviour is still dominated 
by the stable closed loop T{VIIh). As Yj*^ increases to a value of Yj*^ = —.2458, the 
transverse wrapping of the orbits around the closed loop T{VIIh) is increasing. The orbits 
wrap transversely around the closed curve !F{VIIh) so much so, that they catch up with 
themselves and create a new stable attractor that is topologically a torus (which we represent 
as Tiyilh))- The creation of the torus is not especially clear in this view, because the 
interior radius of the torus is extremely small. Note the size of the throat of the torus. 




A tiny region of the parameter space (the loophole) was discovered in the Bianchi type IV model which 
contains no stable plane-wave equilibrium points. We proved the existence of a closed orbit in the loophole 
and provided criteria for its stability. With the use of theory and extensive numerical experimentation, a 
limit cycle was found inside the loophole which acts as an attractor (the Mussel attractor). 

We then studied the Bianchi type VII/j models numerically. In particular, from a local stability analysis 
and extensive numerical analysis we found that at late times 51 — s- and the variables 'freeze' into their 
asymptotic values (in a time scale much shorter than the other dynamical features present) . After this freezing 
has occurred, the system of equations effectively decouples making these models amenable to analytical 
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Figure 5. The figure below displays the dynamical behaviour for values of the parameters 
near the loophole 7 = 1.325 for Bianchi type VII^, non- vacuum tilted models. As T,'^ 
increases from — .2458 to E!J_ = —.2399, the interior radius of the torus increases. The 
closed curve !F{VIIh) lies in the interior of the torus. The future asymptotic behaviour is 
dominated by the attracting torus T{VIIh)- Note how the throat of the torus narrows as 
E^j. increases. As increases from — .2339 to E^ — —.2250, the throat of the torus 
collapses, and the stable attractor becomes the equihbrium point C^{VIIh)- 




investigations. This, in turn, showed the existence closed curves which interestingly had a characteristic 
frequency / = 1/Ti = (1 — (a;))(l + E!^)/ {nVh). From a detailed numerical analysis of the type VII/i models, 
we then showed that there is an open set of parameter space in which solution curves approach a compact 
surface that is topologically a torus. The stability of these late-time attractors seems to be preserved under 
the transition IV 1 — > VII/j. (e.g., the closed curve in the loophole in the type IV models to the torus in the type 
Yllfi models). Finally, we would like to emphasise that the comprehensive numerical integrations presented 
in this paper (in the full Bianchi state space) serve to fully justify all of the claims in particularly the 
'freezing in' process, the existence of the mussel attractor in the loophole and the existence of the torus. 
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Appendix A. Numerical Integrations 
A.l. Bianchi type IV and Vllh. 

A. 1.1. Initial Conditions. The initial conditions are chosen so that the constraint equations are satisfied 
initiaUy. In both the Bianchi IV and Bianchi Vllh cases the initial conditions are chosen so that £_ = 
0.1, Y.12 = 0.01, Sia = 0.01, N = 0.1, A = 0.1 with A 1 for Bianchi type IV, and A = 2/V3 for Bianchi type 
VII/i. The initial values for E4. and E23 are chosen so that the constraint equations are satisfied. Note how 
the initial value of E23 changes as a function of 7. 





^23 




7 = 1 


7 = 


= 1.325 


7 


= 4/3 


-0.7 





48980 





43387 





43349 


-0.6 





55218 





51171 





51144 


-0.5 





59992 





57048 





57028 


-0.4 





63632 





61566 





61552 


-0.3 





66325 





65009 





65000 


-0.2 





68184 





67539 





67534 


-0.1 





69275 





69159 





69158 


0.0 





69635 





69003 





68999 


0.1 





69275 





68016 





68008 


0.2 





68184 





66260 





66247 


0.3 





66325 





63674 





63657 


0.4 





63632 





60154 





60131 


0.5 





59992 





55520 





55490 


0.6 





55218 





49461 





49422 


0.7 





48980 





41357 





41304 



A. 1.2. Type IV: Numerical Integrations. For completeness we present numerical integrations of the dynam- 
ical system for the cosmologically interesting values of 7 equal to 1 (dust), 1.325 (a dust/radiation mixture) 
and 4/3 (radiation). Figures El to El depict some of the results of this integration. The integrations were done 
over sufficiently long time intervals [0, 800], but the plots given here are for time intervals [0, 25]. 

A. 1.3. Type VII^: Numerical Integrations. For completeness we numerically integrate the dynamical system 
for the cosmologically interesting values of 7 equal to 1 {dust), 1.325 (a dust/radiation mixture) and 4/3 
{radiation). Figures 1101 to [TSI depict some of the results of this integration. The integrations were done 
over sufficiently long time intervals [0, 800] but the plots given here are for time intervals [0, 25] or [0, 100] if 
necessary. 

A. 2. Bianchi Type VII/i, !F{VIIh) Throat Attractor. Fixing 7 = 1.325 and choosing different initial 
conditions so that the terminal value of S+ is 'frozen' in between —0.069773 and —0.025, we are able to 
observe some rather interesting dynamics. In this parameter range, we again find that !F{VIIh) is a stable 
attractor. What is perhaps more interesting is the manner in which orbits are attracted to this closed orbit. 
The attractor T{yilh) is shown in Figure ITU 

A. 3. Bianchi Type VII/j, T{VIIh) Torus Attractor. FigurelTBlshows three orbits with initial conditions 
"far away" from the type VII^ plane wave solutions that approach the torus attractor T{VIIh) at late times. 
This figure illustrates the fact that there exists an open set of initial conditions for which the orbits are 
attracted to the Bianchi type VII/i loophole. 
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Figure ITHI depicts the shape and formation of the attractor TiVIIh)- We choose different time intervals 
to illustrate the structure of this attractor. We can see how the curve starts to fill up the torus. Irrational 
curves will eventually completely fill the torus. 
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Figure 6. Type IV: The figure below displays the future dynamical evolution for 7 = 
The asterisk indicates the initial condition. Note how ^ and ^ but 7^ 0. 
this case, the local sink is £{IV). 
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Figure 8. Type IV: The figure below displays the future dynamical evolution for 7 — 1.325. 
The timescale has been increased. The asterisk indicates the initial condition. Note how 
^ and — > 1 or V'^ oscillates, but E_|_ -/^ 0. Further analysis of this "exotic" behavior 
is described in the text. In this case, the local sinks are £{IV), J^{IV), C{IV) and C^{IV). 
If we look closely, one closed loop C{IV) can be seen. 
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Figure 9. Type IV: The figure below displays the future dynamical evolution for 7 = 4/3. 
The asterisk indicates the initial condition. Note how ^ and ^ 1, but E-)_ -/-> 0. In 
this case, the local sinks are £{IV) and C^{IV). 
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Figure 11. Type Yllh- The figure below displays the dynamical behavior for 7 = 1.325. 
Note how ri ^ but E+ 7^ 0. Note the increased timescale in the graph of t vs V'^. Only 
two trajectories were sketched in the first figure to show some of the structure of the future 
asymptotic attractor. Note how both of the trajectories in the first figure spiral into what 
appears to be a point near vi ~ —1, but then one of the trajectories starts to spiral out 
again. This is evidence of the torus attractor T{VIIh)- 
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Figure 12. Type Vllh- The figure below displays the dynamical behavior for 7 = 1.325. 
The timescale has been increased in the graph of t vs V'^. Note how V'^ ^ 1 01 V"^ oscillates. 
Further analysis of this "exotic" behavior takes place in the text. In this case, the local sinks 
are £{VIIh), J-{VIIh), T{VIIh) and C-{VIIh)- This oscillating behaviour indicates the 
existence of the attractor TiVIIh) 
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Figure 13. Type VII^: The figure below displays the dynamical behavior for 7 — 4/3. 
Note how — > and V'^ ^ 1, but S+ -/-> 0. In this case, the local sinks are £{VIIh) and 
C-{VIIh)- Note the increased timescale in the graph of t \s V'^ . Only four trajectories 
were sketched in the first figure to show some of the structure of the future asymptotic 
attractor. 
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Figure 14. Closed loop T{VIIh), T.\ = -0.40. The picture below is magnified to show 
additional details of the structure. 




Figure 15. The figures below show three orbits with inital conditions "far away" from 
the plane waves solutions (initial condition is (S+)o = 0.2) that end up in the type VII^ 
loophole. This illustrates the non-zero measure of the attractor basin of the torus attractor. 
Note how the orbits start to oscillate after the limiting value is reached. 




